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ABSTRACT 


It is possible to determine the phases of structure factors of a non- 
centro-symmetric crystal from an experimental measurement of the 
intensity-differences between pairs of inverse feflections, produced by 
anomalous dispersion. But the determination is ambiguous. An inter- 
esting possibility arises when the group of anomalous scatterers pos- 
sesses a centre of inversion. It is found that in this case, a Fourier 
synthesis calculated by using both the possible values of the phase angle 
will contain, in addition to the peaks of the real structure, an equal num- 
ber of negative peaks at positions related to the former by inversion about 
the inversion-centre of the group of anomalous scatterers. This idea 
has been used to determine the absolute configuration of L-ephedrine 
hydrochloride. 


1. INTRODUCTION 


DuRING recent years the anomalous dispersion effect has been increasingly 
used by many workers to determine the absolute configuration of a number 
of crystals. In all the cases, the structures had to be solved beforehand by 
conventional methods of X-ray crystal structure analysis (which give only 
the structure and not the configuration) and the anomalous dispersion effect 
was used only to determine the configuration of the structure. However, 
the effect can be put to far greater use than mere determination of the con- 
figuration. In a recent paper (Ramachandran and Raman, 1956), it has 
been shown that it is sometimes possible to determine both the structure and 
the absolute configuration of a non-centrosymmetric crystal in a straight- 
forward way by making full use of the effect of anomalous dispersion. It is 
found that the phases a (hkl) of the structure factors F’ (hkl) can be deter- 
mined, but for an ambiguity, from an experimental measurement of the 
intensity differences between pairs of inverse reflections hkl and hkl, pro- 
duced by the imaginary component 4f” of the atomic sacttering factors. 
The nature of the ambiguity as well as methods for its resolution have been 
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discussed and the validity of the method was also tested in the case of the alka- 
loid, L-ephedrine hydrochloride. The phase angles determined by the anomal- 
ous dispersion method were found to agree very well with those calculated 
from known atomic co-ordinates. In this paper, further results which tend 
to confirm the various conclusions previously arrived at are reported. The 
absolute configuration of L-ephedrine hydrochloride, which has been deter- 
mined by an interesting application of the new method, is also presented. 


2. PHASE DETERMINATION FROM ANOMALOUS DISPERSION 


The basic theoretical relations have already been derived in the earlier 
paper and only a short account is given here so as to make this paper self- 
contained. Let the crystal contain one atom or a set of p crystallographically 
equivalent atoms A for which the imaginary component 4f” is significant, 
while for the rest, R, of the structure (N — p atoms) this component is negli- 
gible. Let F,’ and F,” denote the contributions from the real and imaginary 
parts of the scattering factors of atoms A and let Fz denote the contribution 
from the atoms R. These vectors are shown in Figs. 1 (a) and 1 (0) for 
both the reflection Ak/ and its inverse hk/. (In the latter case, the vectors are 
represented by a bar over the symbols.) From these figures it is clear that 


|F)?=|F’|?+ |B,” |? +2] F’| | F,”| cos 0 (1) 














. 1. Diagram showing the different components of the structure factors F and F, 





Anomalous Dispersion Method of Determining Structure 
and 
| F |? =|F’ |? + | Fy” |?—2|F’| |F,"| cos 0 
whence 
4|F|?=|F|?—|F|? 
=4|F’| | F,y”| cos 0 (3) 


where @ is the angle between the vectors F’ and F,”. It is obvious that F’ 
is the structure factor when the anomalous dispersion effect is zero. If a be 
the phase of F’ and a, the phase of F,’ then, 


a=4tn+a,— 0. (4) 


Equations (3) and (4) may be used to determine a from experimentally 
measured 4 | F |? provided the positions of the anomalous scatterers A are 
known. It is not difficult to know the latter, since the anomalous scatterers 
are also usually heavy atoms. 


But there is an ambiguity in the determination. This ambiguity arises 
because @ derived from cos @ [equation (3)] may be positive or negative and 
correspondingly there are two values possible for a, given by 


3. THe ‘‘ ANOMALOUS” SYNTHESIS 


It is interesting to investigate the nature of the Fourier synthesis in which 
both the possible values, a, and a,, obtained by the anomalous dispersion 
method, are used. This synthesis may be called the ‘‘ anomalous ”’ synthesis: 
pan(xyz). It is given by 


1 _+e 7 ' 
Pan (xyz) = aV = = a | F’ (Akl) | {exp ia, (hkl) 


+ exp ia, (hkD} exp — 2mi (hx + ky + Iz) (6) 
and the two values of a are related by 
Gg = 2a, + 7 — ay. (7) 


In general, the synthesis (6) may be shown to contain peaks at the correct 
atomic positions, in addition to a background which may contain spurious 
peaks. An interesting possibility arises in the special case when the group 
of anomalous scatterers possesses a centre of inversion. In this case, 
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the value of a, is 0 or 7 and the two possible values of the phase angles ap 
aand(z—a). If both these values are used, the resulting synthesis is simp) 
the Fourier sine synthesis ps, and is given by 


co +00 +60 
Pan (XYZ) = ps (xyz =2 ztzz | F’ (AkD) | sin a (aki) 


x sin 22 (hx + ky + Iz) (j 


This synthesis may be shown to contain, in addition to the positive peak 
of the true structure, an equal number of negative peaks at positions relate 
to the former by inversion at the inversion-centre of the group of anomaloy 
scatterers. The proofs for both the cases will be given later. 


Thus it is possible to determine the structure and absolute configuration 
from the anomalous synthesis. However it is to be noted that the anomaloy 
synthesis will give no information about those parts of the structure which 
by accident happen to be centrosymmetric about the inversion-centre of the 
anomalous group. In such cases it is necessary to use a certain amount of 
care and judgment and confirm the deductions with other methods like the 
heavy atom method suggested in Section 5. 


This method of determining the structure in its absolute configuration 
seems even more promising than the method of Pepinsky and Okaya (1956) 
using the Patterson sine synthesis (Ps). These authors have investigated the 
nature of the Patterson function obtained from the intensity data alone wha 
some of the atoms are anomalous scatterers. It is found that both the struc 
ture and absolute configuration can be derived from the imaginary pat 
P, (xyz) of the complex Patterson function. However the interpretation of 
the Patterson sine synthesis is much more difficult than that of the Fourier 
sine synthesis. This is obvious from the fact that, whereas the latter consists 
of only the structure and its negative inverse, the P,-function will consist of 
a series of p images of the structure as seen from each one of the p atoms in 
the set A and will contain, in addition, a duplicate of the entire series in th 
form of negative peaks by inversion at the origin of the diagram. 


4. COMPARISON OF THE “ ANOMALOUS” AND THE “ ISOMORPHOUS ” 
SYNTHESIS 


It is interesting to compare the superposed synthesis obtained with the 
anomalous dispersion method with a similar synthesis which may be obtained 
in the parallel case of phase determination by the isomorphous replacement 
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technique. In this case also the phase determination is ambiguous and the 
two solutions a, and a, obtained for a are related by 


a, = 2a, — a,  & 


where a, is the phase of the contribution from the group of replaceable atoms. 
The Fourier synthesis obtained by using both the solutions may be called the 
“isomorphous synthesis,’ pjs (xyz), and is given by 


+00 
pis (xyz) = W z z z | F’ (hkD) | {exp ia, (hkl) 


+ exp ia, (Ak/|)} exp — 2mi (hx + ky + Iz) (10) 


As in the case of the anomalous synthesis, this synthesis also will give peaks 
at the correct atomic positions in addition to a background. In the case 
when the group of replaceable atoms possesses a centre of inversion, the two 
lutions are a and — a and the isomorphous synthesis calculated using both 
the possible values becomes merely the Fourier cosine synthesis (p¢ say). 


F — 


pis (XYZ) = pe (xyz) = + ye + aa x | F’ (hk) | 


x cos a (hkl) cos 2a (hx + ky + Iz). (11) 


It can be shown that this pe function will contain two sets of peaks, both 
positive and the two sets will be related by inversion about the inversion- 
entre of the group of replaceable atoms. Thus pigs will give the structure and 
its positive inverse whereas pg» will give the structure and the negative inverse. 


Obviously a happy example is one where the replaceable atoms are also 
anomalous scatterers. Then the isomorphous method can be used to deter- 
mine values of a and —a and the anomalous dispersion method determines 
values of a and (7 — a). Using both methods in conjunction it is possible to 
fx uniquely the value of a and calculate the electron-density, 


o (xyz =e 4 2 2 - E |F’(kd) | cos « (Aki) 
X cos 2m (hx + ky + Iz) 
22 +co +00 ? 
+ G2 E EF (ik) | sin a (Aki) 


x sin 2m (hx + ky + Iz). 
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Thus the two methods are complementary. The isomorphous replace- 
ment method determines the first synthesis given by equation (12) (i.e., the 
structure and its positive inverse) and the anomalous dispersion method 
determines the second synthesis (i.e., the structure and its negative inverse) 


and the two together add up to give the electron-density p (xyz) (i.e., the 
structure). 


It is also interesting to note that both methods are based on more or less 
the same principle, that phase determination is possible by introducing 
a change in the scattering powers of a few atoms and studying its effect on 
the intensities of the X-ray reflections. In the isomorphous replacement 
method, this change is achieved by isomorphously replacing some of the 
atoms. In the anomalous dispersion method there is no replacement as 
such, but use is made of the fact that when anomalous dispersion takes place, 
the intensities of the inverse reflections hk/ and hk/ are different. The phases 
are determined from the intensity differences between pairs of inverse reflec- 
tions. Thus in this case the isomorphous substances appear to be merely 
the structure and its inverse and the analogy is clearer if it is remembered that 
the reflection hk/ of the structure is the reflection Ak/ of the inverse. 


Between the two methods, preference should be given to the anomalous 
dispersion method because it requires only one crystal, whereas the other 


method demands a series of isomorphous substances, with a minimum of at 
least two. 


5. EXPERIMENTAL VERIFICATION 


The method has been experimentally verified in an investigation of the 
structure of L-ephedrine hydrochloride. The structure of this compound 
has been determined by conventional methods (Phillips, 1954). Our aim 
was to test the new method both as regards the determination of phase and 
the determination of absolute configuration. The crystal belongs to the 
monoclinic space group and the unit cell contains two chlorine atoms. These 
show a measurable, though small, anomalous effect (A f".; = 0-69) even 
for CuKa radiation. The crystals are elongated along the b-axis but the 
hOl reflections are not useful for our purpose. Therefore it was necessary 
to make measurements in the c-zone. The differences in the intensities of 
the pairs of inverse reflections hkO and AkO were measured with a Geiger 
counter spectrometer. The measurements were put on an absolute basis 
by direct comparison with the reported values of Fe. The values of a and 
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(7 — a) were obtained by using an approximate form of expression (3), 
namely, 


sina = 4|F|/2F,’, 


which is valid if F,”<F,’. The phase angles thus derived were found 
to be in good agreement with the correct values calculated from the atomic 
co-ordinates. The measurements were also put on an absolute basis by 
the method of Wilson (1942) which is the usual method adopted when the 
structure is not known. The phase angles derived in this way were not much 
different from the above values. A table showing the phase angles derived 
by the former method and the correct values has been given in the earlier 
paper (Ramachandran and Raman, 1956). 


Using these values and the observed values of | F,| the pan-synthesis 
was calculated (Fig. 2). The phases of some of the reflections could not be 
determined because the contribution from the chlorine atoms to these reflec- 
tions is zero. The corresponding coefficients had to be dropped from the 

















y 

Fic. 2. pay-function projected on (001) using values of phase angles found by the anomalous 
dispersion method. Contours have been drawn at arbitrary intervals of 20. The full lines repre- 
sent positive contours and broken lines represent negative contours. 


synthesis. Further, the phases of a few reflections disagreed by more than 
twenty degrees from the correct phases calculated from known atomic co- 
ordinates. In spite of this, these were included in the synthesis, as this is 
what one would do if the phases were not known to start with. To see what 
would be the effect of both these defects, the synthesis was calculated using 
the correct values of a and (7 — a) taken from Phillips’ paper and this is 
shown in Fig. 3. It is found that the two diagrams agree in their broad 
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features, though the latter is in better agreement with the correct structure. 
In interpreting the diagrams, the important fact should be borne in mind 
that the positive and the negative density distributions are superposed and will 
interfere. This effect will be more serious in a two-dimensional projection 
than in a three-dimensional synthesis. To see the nature of this interference, 
the two molecules are drawn in both the diagrams, using the atomic 
co-ordinates given by Phillips and their negative inverses have also been 
drawn. If now the phases are calculated, taking only the positions of the 
positive peaks then the ambiguity in the solution of a may be resolved. 


4 

















Fic. 3. pen-function projected on (001) using correct values of the phase angles. The con- 
tours have been drawn at arbitrary intervals of 20. The full lines represent positive contours 
and broken lines represent negative contours. 


It may be resolved in another way also (Ramachandran and Raman, 
1956). Of the two possible values of a, the one which is nearer to the phase 
of the chlorine atoms may be chosen for a preliminary Fourier synthesis. 
This procedure is justified for those reflections for which the contribution 
from the heavy chlorine atoms is fairly large. This synthesis is shown in 
Fig. 4, which is seen to be a good approximation to the structure. This use 
of the heavy atom method has been suggested by Peerdemann and Bijvoet 
also (1956). 


6. INDEXING THE REFLECTIONS FOR FINDING ABSOLUTE CONFIGURATION 


The pan-diagram can be used to determine the absolute configuration of 
the compound. The positive peaks of the diagram correspond to the struc- 
ture in its absolute configuration provided correct signs have been assigned 
to the indices hkO. This is easily done by choosing a particular co-ordinate 
system on the crystal and indexing the reflections on the basis of this system. 
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We have conventionally taken this frame of reference to be a right-handed 
system of co-ordinates. For a monoclinic crystal, the senses of the co- 
ordinate axes can be chosen in two equivalent ways to get a right-handed 
wstem, as Shown in Figs. 5(a) and 5(6). The sense of either the c- or a- 
axis can be arbitrarily chosen and then the senses of the other two axes are 
gutomatically fixed by the following two stipulations :— 


(a) The angle between the positive directions of c- and a-axes should 
be obtuse; this fixes the sense of a, once the sense of c has been arbitrarily 
chosen ; 

(6) The sense of 5 is fixed by the fact that the three axes should form a 
right-handed system. 

A 














4 





Fic. 4. Electron-density projected on (001). Contours have been drawn at intervals of 
approximately 1 ¢. A~* except near the chlorine atom where the interval is 2 e. A-*. The one 
tketron cortour is not shown. 
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a Fic. 5. The two equivalent right-handed co-ordinate systems possible for monoclinic 
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In the particular example, the crystal was set with the c-axis vertical. 
The directions of a- and b-axes were determined from the known positions 
at which the corresponding planes (100) and (010) reflect in the spectro- 
meter. The positive direction of c was arbitrarily chosen to be vertically 
upwards. The sense of a can then be fixed from stipulation (a) provided it 
is known whether the chosen direction of a makes an obtuse or an acute angle 
with the chosen c direction. This point was decided by studying the direc- 
tion of shift of the higher layer spots of an oscillation-photograph taken 
about the c-axis. The sense of the remaining axis, namely b-axis, can be 
fixed from stipulation (b). This procedure of assigning correct signs to the 
indices is essential, for a reversal in the sign will lead to the inverse configura- 
tion. Some aspects of the method of assigning the indices correctly in an 
orthorhombic crystal have been discussed by Peerdemann and Bijvoet (1956). 


7. ABSOLUTE CONFIGURATION OF L-EPHEDRINE HYDROCHLORIDE 


The absolute configuration of L-ephedrine hydrochloride has thus been 
determined (Fig. 6). The alkaloid ephedrine CsH;CH (OH) CH (CH,;)NH CH, 
has four isomers, namely, (—)-ephedrine, (+-)-ephedrine, (—)-psi-ephedrine 
and (+)-psi-ephedrine. Their conformations have been studied by stereo- 
chemical methods (Henry, 1949; Close, 1950) and relative configurations have 


Fic. 6. Absolute configuration of L (-++)-ephedrine molecule. 


been assigned to the isomers from their relation, established by chemical 
means, to D-glyceraldehyde and hence to D-tartaric acid whose absolute 
configuration is now knowa (Bijvoet and others, 1951). The present analysis 
gives an independent determination of the absolute configuration of L-ephe- 
drine hydrochloride. The configuration found by us is the one assigned to 
(+)-ephedrine. The earlier X-ray analysis of Phillips did not make use of 
the effect of anomalous dispersion and so it is consistent with this configura- 
tion as well as its mirror-image. Of the two possibilities Phillips chose the 
configuration which was in agreement with the one deduced by a number 
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of chemical workers as the correct configuration. The present analysis con- 
firms this choice. 


The author wishes to express his deep indebtedness to Professor G. N. 


Ramachandran for the guidance and encouragement throughout this investi- 
gation. 
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§1. INTRODUCTION 


THE aim of this paper is to indicate a method of enumeration of positive 
rational numbers. Such an enumeration consists in establishing a (1 — 1) 
correspondence between the elements of the set R+ of positive rationals (or 
a subset of it) and those of the set I* of positive integers. A few methods 
of enumeration of rationals have been indicated in Wilder* (pp. 80-81, 106- 
107). 

In the following, we shall also study some of the properties of the median 
series (defined below) which is akin to the Farey series and also about the 
t-function defined by Hermes.1. We make use of the properties of the 
median series and the 7-function to enumerate the positive rationals > 1. 


Hermes, on the basis of his 7-function, proves certain results which, 
we are bound to think, are motivated by the structure of the median series 
studied in the paper. However Hermes does not obtain all the results 
obtained here and in particular Hermes does not interpret the binary scale 
representation of an integer n in relation to the structure of the median series. 


Hermes’ paper was as a matter of fact suggested by a paper of Stern? 
in which he introduces the series called the Stern’s series 


s.. Gs 
S,.. (r+s,s) 
S,.. (7,2r+s,r + s,r + 2s, s) 


and so on. If r=0,s= 1 in S; then this series becomes identical with the 
numerator of the median series M;. If r=1, s=0O, the series becomes 
identical with the denominator of Mj. From our point of view the import- 
ance of Stern’s series is subuidinate to that of the median series. 


The significance of the partial quotients in the division transformation 
for two numbers a and 3 relatively prime, is also brought out in terms of 
the median approximations defined in Section 5. These approximations 
12 
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are also the approximations of Vahlen® but they differ from those of Vahlen 
in the method of introduction. While Vahlen introduces the approximations 
to a real number p in connection with Farey series and is able to prove the 
important result that the approximations can be arranged in a single series 
ee re so that a;, aj,, are the two consecutive terms of F;, the 
Farey series of i-th order, between which the real number p lies, we however 
define the same approximations in relation to the structure of the median 
series and Vahlen’s series arrangement is no longer true for the median 
approximations defined by us. The relation between the Farey series and 
median series do not seem to be fully known though we can show that the 
Farey series F; is a part of the median series My but not a part of Mjy_,. The 
latter part is evident since the terms 1/i belong to M,; but not to Mj.,. To 
prove that Fyc M,; we can use the following result which can be easily proved: 
‘“* The sum of the partial quotients of a rational number m/n in lowest terms 
is equal to or less than the greater of m and n.” 


Hermes’ 7-function is introduced here in Section 6 specifically by means 
of the order in M, and the recurrence relations on which Hermes bases his 
definition of the 7-function are directly deduced in Section 6, from the order 
properties of the median series. 


The results in Sections 5 and 6 are believed to be new. The auxiliary 
properties of medians, necessary for proving them, have been collected in 
Section 4. 


I am indebted to Dr. R. Vaidyanathaswamy for his guidance and to 
Mr. S. Swetharanyam for his help in preparing this paper in its final form. 


§2. NOTATIONS AND DEFINITIONS 
2.1. The 7-function of Hermes is defined by the recurrence relation 
7, = 1 and t= 7 (n-2”) + 7(2**244-n) 


where 2” < n < 2”+1, 


2.2. The Gaussian bracket is defined by (2.1.1) 
[ }]=1, [a]=—aand 
Ss a <2 Seer (2.2.1) 


=n fa, f,....., Mile, g,..-..G)-+ te &.....,. 9) fiy.....5 80 


2.3. If any positive integer n has the unique alternating binary scale 
expansion 


n= 20 — 20 4 2 + 2%, (2.3.1) 
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(0 < a < a, < a, < .... <a) with odd number of terms (i.e., A even), 
then 


Tn = T (do, 2, — Ag, Ag — A), +2009 Qy — Gy) (2.3.2) 
= (1 + dp, a, — dy, dg — Ay, ....5 A, — Ay-4) (2.3.3) 


[vide Hermes? (satz 1)]. The fact that the alternating binary scale expansion 
is unique for a given integer and that there are an odd number of terms in 
such an expansion is true for n = 1, 2, 3 for, 1 = 2°,2=2!,3=2°-—2 
+ 2%. Let us now assume that this is true for all integers less than N (say) 
where 2-1 < N <2*. Now N can be uniquely written in the form 


N = 2%14 2% 4904 4-90, 


where kK — 1 >4,>4,> .... >Qm=>0. (2.3.4) 


Let now the alternating binary scale expansion of N — 2-1 be 


N— 21 = 2041 904 (odd number of terms) 


Then 
N = 2h-14 Qa gb4 (even number of terms) 
=2.2k-1— Qk-14 Qar1_ 9b 


= 2k — 2h1+4 201204 |... (odd number of terms). 


This is an alternating binary scale representation of N with an odd number 
of terms and this is unique because (2.3.4) is unique. 


2.4. The Median series is defined as follows: 


w=(04) 


m= (oh 


StLizati 
Dares 


mn -GLisigte7 


M. 


and so on. 
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§3. DIVISION TRANSFORMATION AND EUCLIDEAN ALGORITHM 


It is known that given any two positive integers a and 5, an Euclidean 
algorithm can be constructed in the form 


a = ba,+ by (0 <b, < b) 

b = ba, + d, (0 <b, < By) 

by = bya, + bz (0 <b, < by) 

bn—3 = Dn—2 An + On (0 < bn, < bn-2) 

bn-2 = bn an (an > 1) (3.1) 
which ultimately yields the Greatest Common Divisor (hereafter denoted as 
(GCD) of a and 6. The numbers dp, a, de, ...-, @, Which appear in the 


successive stages of the division transformation are called the partial quotients 
in the Euclidean algorithm. It is also known that the rational a/b can be put 
in the form of a regular finite continued fraction (hereafter denoted as CF) 


a l l l 

or a a" eee see ae (3.2) 
The numbers Cp, ¢;, Co, .... are called the partial quotients of the continued 
fraction. It is easily seen that the equations 

a= CG (i = 0, 1, 2,3, ....,) 


are true; that is, the partial quotients obtained in the above two operations 
are respectively the same. 


The number of partial quotients can always be made odd; for if it is 
even, then ay, can be replaced by the two partial quotients (ap_,, 1) so that the 
CF becomes 


a 1 1 1 1 

praterat eee +e + 
In (3.1) by, is the GCD of a and 5. The remainders bo, b,, be, ...., bn 
at each stage of the operation can be eliminated in terms of dp, a, da, ...-5 An 


and by, and then we get 
bn-2 = bn-1 an 
bn-s = Dn [Gn-1, an] 


Dn—s = On-1 [Qn-2, An-1, Qn] 


eevee e eee eee ee er ere e es esereseeesesesers 





= by-y [d2, ds, -- ++» An]. 
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Finally we have 
b = by_, (aj, ag, Az,..--, An] 
and 
A = Dns (Gos A, Ags « - + +5 An) (3.3 
Since bn_, is the GCD of a and 3, it follows that 
the integers represented by the two Gaussian brackets in (3.3) are 


relatively prime. (3. 


The result of the elimination is also the expression of a/b as the regular fini 
continued fraction 





a 1 1 1 

pete + a, tees +e 
ini [4o, 41, Ae, s+ ++9 Qn]. 
[41, de, Ag, .---5 An] 


It follows that if a/b is in its lowest terms 
a= [a0, Ay, Ag, --- +> Qn] 
b = [d,, dg, dg, ..--, An]. (3. 


From the above we see that the positive rational R = a/b can be ma 
to correspond to a sequence of integers dp, d,,...., an (dg > 0, a> Oh 
i= 1, 2, 3, ....,) which are the partial quotients in the CF expansion 
R and these partial quotients completely specify the rational number. 


a/b > 1 (so that dy > 1) is in its lowest terms, then the reverse of the CFi 


equal to another rational number 





R’ _— [an, an-15 seers Ay a] 
[@n-1, An-2) +++ + a, ao] 


= 2 


~ Dns 


where py, is the numerator of the n-th convergent of the original CF. We 
call R’ as the associate of R. It is clear that the relation between a ratio 
number and its associate is symmetrical. 


Therefore we obtain that if a4, >1 (i@=0, 1, 2, ...., a) 


[ao 4, - +++, An] = [Gns Qn-1) +++ +5 Gy, Ag]. (3. 
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With any rational number R = a/b we shall associate the function o (R) = 
the sum of the partial quotients in the CF development of R; i.e., 


o(R)=a,+a,+.... +a. If R and R’ are associates, 


then o(R) = o(R’). The significance of this function is brought out later 
on. 


§4. THE MEDIAN AND THE ELEMENTARY PROPERTIES OF THE MEDIAN SERIES 


We define the median of two positive rationals a/b and c/d (in lowest 
terms) to be (a+ c)/(6+d). Some of the important properties of the 
median are 


4.1.1. The median m,/m, of two unequal positive rationals a/b and 
c/d lies between them and is different from them (a property of ratios). 


4.1.2. The determinants bc — ad, m,b — m,a, m,c — m,d of a/b, m,/m,, 
c/d (taken two by two) are equal. Conversely, if the determinants of the 
above type are equal for three ratios, then each one is the median of the other 
two. 


4.1.3. Ifbc — ad = 1, then the median of a/b and c/d is in lowest terms. 


4.1.4. If m,/m, is the median of a/b and c/d, then m,/m, is the median 
of b/a and d/c. 


For the construction of the median series, we start with the extreme 
irreducible fractions 0/1 and 1/0. We denote M,= (0/1, 1/0). Any 
Mii, i= 90, 1,2, .... is obtained from Mj by itroducing between every 
two successive terms of M,; their median. Thus we get the median series 
mentioned in (2.4). Mj is called the median series of i-th order. The 
following are some of the basic properties of the median series: 


4.2.1. Every term in each of the median series is in lowest terms and 
the terms in a median series are in natural order from left to right and are 
all distinct. Further any two successive terms in a median series have deter- 
minant 1. 


4.2.2. Every term of Mj not in M;_, is the median of the two adjacent 
terms in Mj_,. These terms are irreducible. Every term of Mj; which occurs 
also in M,_, is the median of the two adjacent terms in Mj. These terms are 
however necessarily reducible. 


4.2.3. The median series M; contains 2‘+ 1 terms. Hence the even 
(odd) terms as counted from the left are also the even (odd) terms as counted 
from the right. The terms in Mj;,, that are not in Mj are the even terms and 
these continue as odd terms in succeeding series. 

A2 
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4.2.4. The r-th term from the left in Mj is the reciprocal of the r-th 
term from the right. So the middle term is its own reciprocal and so is 1. 
(It can be seen that the median series of i-th order is got by adding to the 
Farey series of i-th order, its reflection about 1.) 


4.2.5. The middle term of the left half of Mj; is 1/2. The r-th terms 
to the right and the left of 1/2 add up to 1. 


4.2.6. The r-th term from the left (right) in M; occurs as the (2r — 1)-th 
term from the left (right) in Mj;,; hence as (4r — 3)-th term in Mj,, and so 
in general as {2* (r — 1) + 1}-th term in Mj.x. 


4.2.7. So if 2° is the highest power of 2 dividing r — 1, then r-th term 


from right (left) of Mj occurs in Mj-, as the 2 a-th term from the right (left) 
where 


r—1l 
a= eq + t- 


Being an even term, this term occurs for the first time in Mj_¢. 


§5. THE MEDIAN SERIES, (VAHLEN’S) APPROXIMATIONS AND THE 
ENUMERATION OF POSITIVE RATIONAL NUMBERS 


The following is an interpretation based on the idea of approximations 
to a rational number in terms of the partial quotients ap, a, dg, ...., dn of 
the CF expansion of the rational number and this brings out the relation 
between the approximations and the median series and also the fact that the 


rational number p occurs for the first time (and so as an even term) in M,i,). 


It is clear that M, has no term betewen pand 1/0. It can be seen that the 
largest value of i such that there is no term between p and 1/0 in Mj is a. 
The penultimate elements 1, 2, 3, ...., @) in M,, Mg, ...., Ma, respectively, 
may therefore be called the first a, median approximations to p.* 


Again the largest value of i such that in Mg,,;4 there is no term between 
p and a,/1 is evidently a,. The numbers ay + 1/1, ay) + 1/2, ...., ag + I/a, 
Of Ma,i1, Mae, ----> Maia, respectively, in the interval (p, 1/0) give the 
second set of median approximations to p.f Proceeding thus we have 


5.1. ais the maximum value of ¢ such that there is no term in 


Ma,+0,+0;+- ee oe between p and ay + l/a, + 1/a, +....+ 1 /aj-, 
@=Z, 3, .... Rp 


* The successive terms in M, between which p lies are (i/1, 1/0), i= 1, 2, 3, ...., a 
+ The successive terms in Moos, between which p lies are (ao/1, a) +1/i), i=1, 2, 3, ...., ay 








—_— 
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5.2. The ith set of approximations to p are: d+ 1/a,+ l/a,+.... 
+ Ifa, + Aft (t=1, 2, ...., @iy;s i= 2, 3, oesg tt ER 


5.3. Thus the total number of approximations to p is ad) + a,+.... 
+ Gn =«(p). This shows that p occurs for the first time in M,,,). It is 
to be noted that all the convergents of the CF occur among these approxima- 
tions. 


We now proceed to find the exact order in which the convergents dp, 
My + 1/ay, 4g + 1/a, + 1/ag,.... occur in the respective median series Mg,, 
Mg,ia,. ----- From the above reasoning relating to the approximations 
it is clear that 


5.4. d, and 1/0 are consecutive in Mg, and p lies between them; 


5.5. d,and a)+ 1/a, are consecutive in Ma,;a, and p lies between them ; 
and so on. We note that a) occurs in Mg, as the 2%-th term and hence occurs 
in Ma,+a, as {2% (2% — 1) + 1}-th term. Since a+ I/a, is the term next 
to dy in Mg,sa, its order is 2% (2% — 1)+ 2. Hence a) + 1/a, occurs in 
Ma,+a,+a, aS the {2% (2%+4 — 2% + 1) + 1}-th term. Since a) + 1/a, + 1/az 
is the term preceding a) + 1/a, in Ma,+a,+, its order is 2%+% —20+4+ Jas, 
Hence in general 


5.6. the i-th convergent a,+ l/a,+ 1/a,+....+1/a occurs in 
Ma,ia,+----+a, a8 the p-th term where p = 2%+@+----+0;_Qas+... +061 
seve (12 + {1 + (—1)}. 


Proof.—Let this be true for the i-th convergent. We shall prove that this 
is then true when i is replaced by i+ 1. The term (—1)** indicates that if i 
is old, the order of (i + 1)-th convergent in Mg,;...... sai¢, JS ON€ more and 
if i is even one less, than the order in Mg,,....ais. Of the i-th convergent. 
From (5.6) and (4.2.6) the i-th convergent has in Mg,;...-+ai+:, the order 


2941 (p — 1) + 1=2%,.. 1041 |.) 4 (— 1) 2041 + 1, 
Hence the (i + 1)-th convergent has in Mg,,...-.+ai+. the order equal to 
ye ae (—1)i Qaita 4+ {1 en (—1)+%}, 


which is the same as (5.6) with (i + 1) instead of i. Since the result has been 
verified for i= 1, 2 it holds generally. 


The above is also true for the n-th convergent, that is p. For thé sake 
of convenience let us take n to be even so that the CF has an odd number of 
terms. We then have the 
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THEOREM |. A rational number p occurs for the first time as q-th term 
in Map); if p has partial quotients ao, a,, ...., an (n even) then, 


Qe Ft 9 A 8 J A, + 2a 
~ 11 1 
PPrRtS Tet + HZ. 


It can be seen that since ay > 0, g is even as ought to be the case since 
p occurs in M,,,, for the first time. 


COROLLARY.—The number p occurs as the 2r-th term of Mj», where 


o=a(p)=—at+a,+.... +a (n even) 
ret — Qe + 21, 


Note.—If p < 1, then a, =0 and in this case the first two terms of r 
cancel each other. 


From the corollary we see that we can associate with the partial quotients 
of p, the integer r= 27-1 — 27-1-@ + .... +21, ¢ = a(p) and then we 
are in a position to say that p occurs as the 2r-th term in Mg, . Conversely 
if we pick out any term in M, which occurs for the first time, its order is an 
even number 2r and we may show as follows that any such r can be repre- 
sented as an alternating binary scale representation with odd number of 
terms starting with 27-1. 


For, we know that any integer r can be uniquely expressed as an alter- 

nating binary scale representation with odd number of terms in the form 
r= 2 — 2 4+ 1... + 2 (n even) 
G>aq> a> .... >a =D. 

If r > 2+? (so that p> 1) its alternating binary scale representation will 
begin with 2*1 and we can, by the corollary, obtain its partial quotients and 
in particular the first partial quotient of the corresponding 2r-th term in Mj 
will be greater than zero. If r < 2*-* (so that p < 1), the alternating binary 


scale representation will begin with 2™ (m <i — 2); we shall then add to 
this representation the initial term so that 


r=21—Qeii gm ...,, 


The partial quotients are now given in order (0, i— 1 — m,....). 


5.7. Thus we are not only able to give the precise order 2r of the 
number p and the median series M, in which it occurs for the first time but 
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conversely, given any term, say 2r-th in an arbitrary Mj, we are able to obtain 
its partial quotients. 


Consequently, the corollary establishes a (1 — 1) correspondence bet- 
ween a given rational number p > | and an integer r in the form of an alter- 
nating binary scale representation with odd number of terms, making it 
possible to enumerate the rational numbers greater than 1. 


We next deduce from the corollary, the following theorem which shall 
be used in the next section. 


THEOREM 2. The sequence of numbers which are the R-th terms from the 
right in all the median series, form an increasing arithmetic progression with 
common difference 1. 


Proof 
Case 1.—Let R be an even number 2r. Let p be the 2r-th term from the 
right in Mj and let it have the partial quotients a), a,, ...., a, (neven). Since 


2r-th term from the right is the (2 — 2,r + 2)-th term from the left, we have 
by Theorem 1, 2(241—r+ 1) = 2'— 2+ 4 .... + 29 (7 even) so that 


r= 2-1 _ 7-1-0 + ccs Sits | (5.8) 
This is therefore the unique alternating binary scale representation for r. 
If p’ be the 2r-th term from the right in Mj,, similarly we have, 
r= Ji+k-1-bo alte Jit k—1—be-bi + age Qbm-1 + 1 


where bo, b,, ...., Dm are the partial quotients in the CF expansion of p’. 
Comparing this with (5.8) we have 


bo =a tk; H=a(i>90); m=n. 


Thus the terms in M; and Mj, differ only by k, the difference in the orders 
of the median series. 


Case 2.—Let now R be an odd number = 2°j + 1 (j odd). Then by 
(4.2.7) it follows that R-th term of Mj occurs as (j + 1)-th term from the 
right in Mj-¢ and this is even. Similarly R-th term from the right in Mjix 
occurs as (j + 1)-th term from the right in Mj,z-¢. Since j + 1 is even, the 
arguments of case (i) hold here and we have the theorem. 


From (4.2.4) we now have 


5.9. If the R-th term from the left of Mj is a/b, then the R-th term from 
the left of Mix is (a/ka + 5). 





M. S. SRINIVASAN 


5.10. Hence, the numerators of the R-th terms from the left in My, 
Mi, --.-, (R < 2**) are the same; that is, the numerator is independent 
of i. Also, the denominators form an arithmetic progression with common 
difference equal to the corresponding numerator. 


§6. THE RELATION BETWEEN THE TERMS OF THE MEDIAN 
SERIES AND THE 7-FUNCTION 


The corollary under Theorem | and (5.7) show that the rational numbers 
are completely enumerated by the sequence of the median series. If any 
rational number p occurs for the first time in M, as 2r-th term, then we call o, 
r as the median co-ordinates of p. From what has preceded, the partial 
quotients of p can be found by the expression of r in the alternating binary 
scale. p itself can be expressed by Gaussian brackets in terms of the partial 
quotients. 


But we can adopt a short method by expressing the numerator and the 
denominator of p directly as a function of the median co-ordinates. We 
may define the numerator of the 2r-th term from the left in Mj to be (r, i). 
But from (5.10) it follows that u(r, i) is independent of i; that is, u(r, i) 
=(r) say. This is in fact the r-function of Hermes defined in Section 2. 
The relation between these two functions is given by w(r) = 77. Also the 
denominator of the 2r-th term in Mj is the numerator of its reciprocal, that 
is, the (24 -— 2r — + 2)-th term. Therefore the denominator of the 2r-th 
term 1S 


pe (22 — r+ 1) = Tyger 
Thus 


2r-th term of M; = —— (r <2+-», 
(2 T+1) 


and this is in lowest terms. 
We now proceed to show a method of calculating ry. 


We have seen that given a positive integer r which has the unique alter- 
nating binary scale representation 


p= Qi — 2-1 4 Qeddeds _ | +. 245-1 (n even) 
(a) > 9, dn = 1), 
it follows from (2.3.1) and (2.3.2) that 
Tr = [Gn, Anz, --- +5 Ag] 
= [do, a, ... +, An] by (3.6) 


= [a2, ag, seers An) if dg = 0. 
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From the alternating binary scale representations of the integers 1, 2, 3, 4, 5, 


VizZ.5 
1=2°; 2=2!; 3 =2?— 2!+ 29; 4= 23; 5 = 23 — 2? +2°, 


we get 
~1=1, 7e=2, —=3, m4=3, —%=4. 
The recurrence relation for t, can be obtained from (5.9). For, from (6.1) 
2r-th term from the left in 


Tr 


M; = 
T(9t—1~r+1) 


(r < 2-1) 


Hence 


2r-th term in 





Tr TP Toya) 
As these have been remarked to be in lowest terms we get, 
Taira) = Tr + Toy) (r <2") 


Replacing i by (i+ 1) and the r by (2*— a+ 1) we get 


Tia) = Ta + T2!-at1) (@< 2') (6.2) 


This is precisely the recurrence relation of Hermes! (p. 372). 


We can now calculate any term (say k-th from the left) of any My. Let 
the term be denoted by T (f, i). 


Case 1.—Let k be even = 2r. 


Then 


T (k, i) = —* 


T (28 4-441) 


by (6.1) 


a. 
T42'-k+2) 
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Case 2.—When k is odd, k + 1 and k — 1 are both even. Since k-th 
term of Mj is the median of the adjacent terms, we have 





T (k, i) = —THe-n 1 Tae 6.3 
(Ki) Ty 2'-Kk+8) + 7h 2!-k+1) (6.3) 


Or alternately, when k is odd, write k = 2°j7+ 1 (jodd). Then this term 


occurs for the first time as (j + 1)-th term in Mj_¢, i.e., as the (j +1)/2-th even 
term. 


Hence T(k,i)= T(j+1,i—c). Hence by case (1) we have 


T (k, i) = Gee where k = 2°j+1. (6.4) 


TH 2 *—J-+1) 
It can be verified that (6.3) and (6.4) are the same, using 7-function proper- 
ties. 
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1. INTRODUCTION 


EXPERIMENTAL studies of the elastic properties of single crystals are of great 
significance from the standpoint of solid state physics. Such studies, when 
extended to higher temperatures, are likely to throw light on the nature of 
interatomic and interionic forces in crystals. Hitherto, though it has been 
possible to obtain data for quite a large number of single crystals at room 
temperatures, extensive work on their temperature variation was not done. 
Hence it was thought worthwhile to study the temperature dependence of 
the elastic moduli of some cubic crystals. 


Eight cubic crystals, elastic constants of which have already been studied 
at room temperature by Bhagavantam and his co-workers, have been chosen 
for the investigation in the temperature range 30-300° C. Some of the more 
important results have been published in a recent review article in these 
Proceedings (Bhagavantam, 1955). 


2. EXPERIMENTAL DETAILS 


The composite piezoelectric oscillator method, originally developed by 
Balamuth (1934) and Rose (1936) has been employed in the present investi- 
gation. In this technique, the specimen crystal in the form ofa bar is 
cemented, with a suitable adhesive, to a quartz rod of identical cross-section. 
The fundamental resonant frequencies of the quartz and the composite sys- 
tem are determined separately and the natural frequency of the specimen 
is computed therefrom. The method of finding out the frequencies makes 
use of the fact that at resonance the impedance of the system undergoes 
a marked reduction. With a knowledge of the resonant frequency, length 
and density of the specimen, one elastic constant or one combination of 
elastic constants, can be calculated. 


Quartz rods of 0-30 cm. square cross-section are used to excite the longi- 
tudinal mode. Torsional vibrations are excited by a Y-cut quartz rod of 
circular cross-section and 0-30 to 0:40 cm. in diameter, with the electrode 
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arrangement as suggested by Giebe and Scheibe (1928). The composite oscil- 
lator with its holder is placed at the centre of an electric furnace and the 
temperature is controlled by varying the heating current. 


The lengths of the specimens are measured at room temperature on 
comparator. The densities are determined by the hydrostatic method. 
Densities and lengths at higher temperatures are evaluated using the data on 
thermal expansion, available in the literature. 


Young’s modulus, Y, and the rigidity modulus, n, along (100), (110) 
and (111) directions are directly determined and the elastic moduli are 
evaluated using the formule :— 
1/Y (100) = S113 1/” (100) = Sa4 
L/Y (roy = $ (Sra + Sa2 + $544)3 1/M G0) = Si — S12 + 4 Sa 
1/A ayy = 4/3 (Su — S42) + FS aa. 

The elastic constants, c’s are computed using the well-known formule. 


The estimated uncertainty in s,, and sy, is one per cent. A small per- 
centage error in measurement is very greatly magnified in the calculation of 
the difference quantity s,,. The uncertainty in s,. varies from 2 to 3 per 
cent. over the temperature range. With this method, changes in moduli can 
be more accurately measured than their absolute values. 


a 








RESULTS AND DISCUSSION 



















(a) Alkali halides—Among the cubic crystals, the alkali halides consti- 
tute an important group whose physical properties have been investigated 
extensively. Direct theoretical evaluation of the elastic constants has also 
met with some success in this case. 


The effect of temperature on the elastic moduli of KCl and KBr has 
been studied up to room temperature by Durand (1936) and Galt (1948) 
respectively. In the present investigation these studies have been extended 


to 300° C. and LiF, which belongs to the same type of structure, has been 
freshly studied. 


Ammonium chloride is an interesting case among alkali halides. It has 
a body-centred CsCl structure at room temperature and changes over to the 
NaCl type at higher temperatures. Thermo-2lastic behaviour of ammonium 


chloride is studied with a view to compare it with that cf the NaCl type 
halides. 


In Table I, two constants characterising the thermal behaviour of crystals 
in longitudinal and torsional modes are given. The constant represents the 
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variation of the elastic modulus with temperature expressed as a fraction 
of its value at room temperature. 


The data on NaCl (Subrahmanyam, 1954) arealso included to facilitate 
a comparison of alkali halides. The s’s and the c’s are respectively expressed 
in units of 10-'%cm.?/dyne and 10" dynes/cm.? 
TABLE I 
Thermo-elastic constants of alkali halides 








Substance ~ q ow x 105, x —_ x 105 
NaCl $a 155 35 
KCl ‘a 93 30 
KBr ‘i 113 22 
LiF * 92 27 
NH,Cl ba 40 172 





Experimental results reveal the following points regarding the thermo- 
elastic behaviour of alkali halides. 


1. Rate of variation of s,, with temperature is greater than that of sg, 
in the alkali halides of NaCl type. 


2. Contrary is the case with ammonium chloride as the variation of S44 
is greater than that of s,, in this case. Ammonium chloride belongs to the 
CsCl type structure. 


3. Sy, decreases with increasing temperature in the NaCl type halides. 
It increases in the case of ammonium chloride. 


4. Elastic constants ¢,,, Cy, and C4, all decrease with increasing tempera- 
ture in all cases. 


5. In the case of thermo-elastic constants for the shear mode, there is 
a striking difference between the alkali halides of NaCl and CsCl types. The 
value of the constant for NaCl type halides is about 30 while its value in the 
case of ammonium chloride is 172. In other words, the rate at which Sq, 
increases or c,, decreases with temperature in the case of ammonium chloride 
Ss very much greater than that in the case of NaCl type halides. This can be 
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explained as due to the low melting point of ammonium chloride and cor- 
respondingly high values for the NaCl type halides. As the melting point 
is approached the lattice loses its resistance to shearing more rapidly 
than at temperatures far removed from the melting point. 


(b) Other crystals——The results of measurements made on Iron pyrites, 
Fluorite and Galena are presented in Table II. 








TABLE II 
Temp. 1 d(5y). ,p5 1 4 (Sas) 5 ae 
Substance range °C. 3, @t * 10 rm , x10 Variation of sy, 
FeS, “ 30-300 17 23 +0-2 to +0-35 
CaF, aa 30-300 35 39 —1-7 to —1-80 
PbS a 30-200 45 65 —2-0 to —2-30 





Iron pyrites crystallises in the cubic system and belongs to the Tp class. 
Its elastic constants have been determined by Voigt (1928) using static 
methods. He reported a positive value for s,.. Later on, dynamic methods 
have been applied to study its elastic behaviour by Bhagavantam and Bhima- 
senachar (1944) and Doraiswami (1947) amongst others. In this work, 
crystal rods are cut from a well formed two-inch cube. The results clearly 
show a positive value for s,.. This positive sign persists throughout the 
temperature range (30-300° C.). Pyrites is unique among crystals to exhibit 
such a peculiar elastic behaviour. s,, and s4, vary almost linearly with 
temperature and sj, practically remains constant in all the three cases. 


The results for magnetite are presented in Fig. 1. Magnetite crystal'ises 
in the holohedral class of the cubic system. It possesses a curie point at 
575° C. and the symmetry of the crystal changes from cubic to orthorhombic 
at 119°C. Elastic constants of magnetite vary with temperature in an un- 
usual manner. Elastic constants of solids generally decrease with increasing 
temperature provided there is no change of state. However, as can be seen 
from Fig. 1, Young’s modulus along (100) direction decreases as usual while 
the rigidity in the same direction increases. The increase is definite, though 
small. Along (111) direction the behaviour is exaclty the reverse. Rigidity 
decreases with increasing temperature while the Young’s modulus increases. 
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The composite piezoelectric oscillator method has been employed to 
study the temperature dependence of the elastic moduli of single crystals of 
KCl, KBr, LiF, NH,Cl, FeS,, Fe,O,, CaF, and PbS in the temperature 


range 30-300° C. 
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Fic. 1. Magnetite. 
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ABSTRACT 


Expressions have been obtained for the Geometrical Sensitivity 
as a function of the zenith angle of a cosmic ray telescope comprising 
of counter trays of rectangular dimensions. The Radiation Sensitivity 
and Cumulative Sensitivity have also been calculated, assuming a zenith 
angle attenuation of the form Ig = Ig.cos?@ for cosmic ray intensity. 


I. INTRODUCTION 


MEASUREMENTS of cosmic ray intensity in the vertical direction are mostly 
carried out by Geiger counter telescopes which comprise of a series of counter 
trays of rectangular dimensions situated one above the other. All counters 
in the same tray are connected in parallel and the pulses from the various 
trays are fed to a coincidence unit which gives a sizable pulse in the output 
only if pulses are fed to all its inputs within a small time interval, termed as 
the “‘ resolving time ’’ of the coincidence unit. The coincidence output rate 
corresponds to single charged particles passing through all the counter trays 
except for (1) a contribution due to showers where different trays could be 
triggered by different particles all belonging to the same shower and hence 
arriving almost simultaneously and (2) a contribution due to chance coin- 
cidences where different trays could be triggered by different particles appa- 
rently unrelated to each other. For coincidence units with resolving times 
of the order of a few microseconds as is usually the case and for more than 
two counter trays in coincidence, the contribution of chance coincidences is 
negligible. The contribution of side showers depends on the material in 
the near environment of the telescope and on the separation of the trays. 
Under optimum conditions this does not exceed about 7% for triple coinci- 
dence telescopes, as investigated by Greisen and Nereson.! Thus the vast 
majority of counts registered by such a counter telescope are due to cosmic 
radiation incident in the cone within which a single particle can cross all the 
trays. 
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A vertical meson telescope is characterised by (a) its dimensions which 
may be given as length, breadth and separation of the end trays or, alter- 
natively, the semi-angles subtended by the telescope in two mutually per- 
pendicular planes and the separation of the end trays, (b) the amount of 
absorber used in between the top and bottom trays. Once these character- 
istics are known it is possible to evaluate the expected counting rate of the 
telescope from known values of the vertical flux of cosmic ray intensity at 
the place of observation and the zenith angle dependence of the same. 


A vertical meson telescope offers a maximum sensitive area for particles 
coming in the vertical direction. A simple geometrical consideration shows, 
however, that the solid angle available for particles coming in inclined direc- 
tions is greater than the one available for those coming in a vertical direc- 
tion. Taking into account both these factors one can calculate the Geometri- 
cal Sensitivity G.S. (6) of any telescope arrangement for different values of 
the inclination 9 which the incoming cosmic ray trajectories make with the 
vertical. Parsons? has obtained an expression for the same in case of a 
meson telescope of cubical symmetry which is recommended as standard 
equipment for the current International Geophysical Year. His method 
can, however, be extended to counter telescopes of any dimensions. We 
have attempted here to get a general expression for telescopes of any given 
dimensions. 


II. GEOMETRICAL SENSITIVITY OF A VERTICAL MESON TELESCOPE 


Consider a geometrical arrangement (Fig. 1) in which the top and bottom 
trays are represented by ABCD and A’B’C’D’ respectively. Let the length 
AB = A’B’ = CD = C'D’ = / and breadth BC = AD = B'C’ = A'D’ = d. 
The separation between the two trays is AA’ = BB’ = CC’ = DD’ =a. 

If the ratios //a and d/a are denoted by 4 and 3 respectively, the semi- 
angles of the telescope in the two vertical planes one along the length and 
the other perpendicular to it are tan-! (4) and tan~ (8) respectively. The 
most inclined direction the telescope can record has a zenith angle @ = tan-! 
VP? + d?/a = tan“ »/4? + §2 and corresponds to rays parallel to the diago- 
nals AC’, BD’, CA’ or DB’. 


Consider cosmic rays incident along a zenith angle @ and an azimuth a, 
where a = 0 for a direction parallel to the breadth ‘d’ of the counter trays. 
The beam will then be incident on the rectangular area OXD’Y of the lower 
tray where 


OX = (/—a.tan 9.sin a) and 
OY = (d—a.tan @.cos a) (1) 
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Fic, 1 
The effective cross-sectiona] area of the telescope normal to the beam 
would be 
A = cos 6.(/—a.tan @.sin a).(d—a.tan 9.cos a). (2) 
Consider the shaded area (Fig. 2) normal to the direction of incidence 
of a beam which has a zenith angle ranging from @ to (@ + d@) and azimuth 
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from a to (a+ da). The solid angle subtended by this beam at the centre 
of the telescope is given by 


R.sin 0-da).Rdé 
dw wl — 


= sin 6.d0.da (3) 





The relative Geometrical Sensitivity of the telescope for directions con- 
fined to zenith angles between 6 and (6 + dé@) and azimuth between a and 
(a + da) is given by 


G.S. (0).d0.da = A.sin 0.d6-da (4) 


The total relative Geometrical Sensitivity of the telescope corresponding 
to a zenith angle @ is 


G.S. (8) = f'A.sin 0.da 


= sin 9.cos 6 [I d-(a, — a,) — a-tan @.1.(sin a, 


— sin a,) + a.tan 0.d.(cos a, — Cos a) 


- (a?/4).tan?0.(cos 2a, — cos 2a,) (5) 
Putting (d/a) = 8 and (I/a) = 4, 
ee — ee ee 





— 4.tan 6.(sin a, — sin a,) + 5.tan 0.(cos a, 
— COS a,) — (tan? 0/4).(cos 2a, — cos 2a,)] (6) 


where a, to a, is the range of azimuth for which the effective cross-sectional 
area A has positive values. 


Assuming that the breadth ‘d’ is less than the length ‘/’, the zenith 
angle @ can be grouped into three ranges for each of which a set of values 
of a, and a, would be effective. 


Case 1.—0 < 9 <tan(8). For this range, the azimuth values range 
from 0 to 27. Since, however, all the quadrants are symmetrical, it is suffici- 
ent to integrate from 0 to 7/2. Thus a, = 0 and a, = 7/2. 


Case 2.—tan- (8) < 9<tan-'!(4). For values of 6 in this range 
(see Fig. 3) the particle would miss the telescope for azimuth a =0. The 
A3 
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cross-sectional area A will first become positive for values of azimuth given 
by 


“ce a en = 
COS a, = H = oa” (d/a):cot 6 = 8.cot 8 


Thus the lower limit is a, = cos-!(8.cot 4). 


The upper limit for integra- 
tion is still 7/2 as in case 1. 
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Case 3.—tan-! (4) < 6 < tan 1/4? + 8%. For this region (see Fig. 4) 
the particles will have a positive cross-sectional area for values of a between 
a, and a,, where 


a, = cos! (8.cot 4) and a, = 2/2 — cos! (4.cot 9). 



































Zenith Angle Response of a Vertical Meson Telescope 35 


For isotropic radiation incident upon the telescope, G.S. (4) is directly 
proportional to the number of particles coming at a zenith @. If, however, 
there is a zenith angle dependence of the cosmic radiation of the form 


I, = Iy-cos* 0 (7) 
the Radiation Sensitivity R.S. (6) is given by 
RS. (6) = cos* 6.G.S. (9). (8) 


III. RESULTS 


In Figs. 5, 6 and 7, the relative Geometrical Sensitivity and Radiation 
Sensitivity for various values of the semi-angles tan-!(8) and tan-!(4) of 
a telescope are plotted. The value chosen for A is 2. 
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Fic. 5. Zenith angle dependence of Geometrical Sensitivity and Radiation Sensitivity for 
a telescope having a semi-angle of 20° in one plane and 5°, 10° and 20° in a perpendicular plane. 
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a telescope having a semi-angle of 45° in one plane and 5°, 10°, 25° and 45° in a perpendicular 
plane. 
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Fic. 7. Zenith angle dependence of Geometrical Sensitivity and Radiation Sensitivity fora 
telescope having a semi-angle of 60° in one plane and 5°, 20°, 40° and 60° in a perpendicular 
plane. 


lar The total counting rate of a telescope is given by 
tan (A? + 3 
N = 4], - J R.S. (9-40 (9) 
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The percentage contribution to the total rate of particles confined to 
the zenith angles between 0 and any value @, is given by the Cumulative Sensi- 
tivity C.S. as 


n-1 4/ A?+82 


C.S. = [100. IRS. (#).d6] = [f RS. (8).d0] (10) 


The Cumulative Sensitivity for various zenith angles and for various 
values of tan“! (8) and tan-!(4) is shown in Figs. 8, 9 and 10. 
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Fic. 8. Percentage Cumulative Sensitivity of a telescope having a semi-angle of 20° in one 
plane and 5°, 10° and 20° in a perpendicular plane. 


A striking feature revealed by the plots of Cumulative Sensitivity is 
that the bulk of the radiation is confined to comparatively small zenith angles 
in spite of large opening of the telescope. Thus, for a telescope having semi- 
angles as large as 60° in both the East-West and North-South planes, about 
70% of the radiation is confined to zenith angles less than about 35°. The 
bias towards smaller zenith angles is increased still further if one of the semi- 
angles of the telescope is small. For example, a telescope having semi-angles 
20° and 60°, counts 70% ofits counting rate within zenith angles of 0 to 25°. 


In designing a telescope for a particular investigation, the research worker 
is interested in the directional sensitivity of his instrument. The present 
work enables the investigator to decide the shape of a counter tray which 
would satisfy requirements of directional sensitivity in the most appropriat¢ 
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way. If, for example, the requirements specify an 80% response due to 
particles within 10° inclination with zenith, irrespective of azimuth of direc- 
tion of arrival, a counter tray of square geometry obviously provides the 
optimum solution. If, however, we are interested in a narrow angle of 5° 
in the E.-W. plane, a rectangular tray of dimensions | : 4 in the E.-W. and N.-S. 
planes would still provide an 80% response restricted to 10° with zenith. 


The authors are grateful to Prof. V. A. Sarabhai for helpful discussions 
and to the Atomic Energy Commission of India for financial assistance. 
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From Mahonia nepalensis D.C., Chatterjee! reported the isolation of an 
yellow base umbellatine, C,,H,,O,N, m.p. 205-07° (decomp.) and a red base 
neprotine, C,s,H,,0,N. Earlier Chatterjee had reported the isolation of 
umbellatine from Berberis umbellata? and Berberis insignis.* On the basis 
of degradation work umbellatine was assigned the structure I* and neprotine 
the structure II.§ 
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Both these structural formule have unusual features and a re-examina- 
tion of the alkaloids of Mahonia napalensis appeared desirable. Following 
the procedure of Chatterjee two bases were isolated by us, a yellow base, 
m.p. 205-07° (decomp.) and a red base, m.p. 206° (decomp.). 


The yellow base was characterised by preparation of a crystalline sul- 
phate and a picrate whose melting points were in good agreement with those 
of the corresponding derivatives of Chatterjee’s umbellatine. Analysis of 
the yellow base after drying at 110° in vacuo gave values in agreement with 
the molecular formula Cy>H,,0;N, 3 H,O. Reduction of the ‘ hydrochloride ’ 
of the yellow base in alcohol solution in presence of Adams catalyst yielded 
41 
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a hydrochloride, m.p. 215-17°, from which a colourless base, m.p. 172°, 
analysing for the formula C.,H.,0,N was obtained. The identity of the last 
compound with d/-tetrahydroberberine was established by mixed melting point 
determination with an authentic sample. Reduction of the hydrochloride 
of the yellow base with lithium aluminuim hydride yielded a base, m.p. 164°, 
identical with dihydroberberine (m.p. and mixed m.p.). The yellow base 
gave an acetone compound, m.p. 168°, identical with berberine acetone com- 
pound (m.p. and mixed m.p.). A sample of berberine regenerated from ber- 
berine sulphate (B.D.H.) after one crystallisation from alcohol yielded crystals, 
m.p. 145° (decomp.). A second crystallisation from alcohol gave crystals, 
m.p. 205-07°.* The latter gave no depression in melting point with the 
yellow base from Mahonia nepalensis. It is therefore proved that umbel- 
latine is identical with berberine. 


It is curious that the m.p. 205-07° has not so far been recorded in litera- 
ture for berberine. Katti and Shintre*® isolated from Coscinium fenestratum 
Coleb an alkaloid, m.p. 205°, whose identity was not established. Pillai 
and Varrier’ in a reinvestigation of this plant obtained an alkaloid whose 
identity with berberine was established by conversion to the acetone com- 
pound from which the form of berberine, m.p. 145°, was regenerated. It is 
also strange that melting or decomposition points, have not been recorded 
for many of the salts of berberine. It is true that berberine and its salts 
undergo considerable colour change before melting, but it was observed in 
the present work that the hydrochloride, the sulphate and the picrate all melt 
with decomposition after a considerable preliminary darkening of colour. 
The melting points of berberine and its derivatives and the corresponding 
compounds from the yellow base of the Mahonia nepalensis are recorded 
in Table I. 


Chatterjee has reported that the red base obtained from Mahonia_nepalen- 
sis charred without melting. From the present work it was found that the 
red base mclted at 206-08° (decomp.) and yielded a picrate, m.p. 210-11° 
(cf. Chatterjee 208-10°). The red base analysed for the formula C.,H.,0;N, 
3 H,O. Reduction of the ‘ hydrochloride’ of the red base proceeded with 
absorption of two moles of hydrogen, yielding a tetrahydro base C,,H,,O,N, 
m.p. 206-07°, which after methylation with diazomethane, yielded a com- 
pound m.p. 146°, identical with d/-tetrahydropalmatine. Ethylation of the 


* Prof. K. Venkataraman (private communication) obtained a base, m.p. 205-06°, from 
Coscinium fenestratum, which gave an acetone compound identical with that frcm berberine. He 
also noted that berberine from commercial berberine sulphate melted at 145° after one crystallisa- 
tion and at 205-207° after a second crystallisation frcm alcchol., 
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tetrahydrobase with diazoethane yielded an ethyl ether, m.p. 124°. A sample 
of -corypalmine was ethylated, oxidised with iodine and reduced with zinc 
and hydrochloric acid to yield dl-O-ethyltetrahydrojatrorrhizine, m.p. 124°, 


identical with the ethyltetrahydroneprotine above. 


It is thus established 


that neprotine is identical with the known alkaloid jatrorrhizine. The rele- 
vant melting point data is presented in Table II. 























TABLE I 
The yellow base Umbellatine Berberine 
(present work) (Chatterjee’s work) m.p. 
m.p. m.p. 
Base .. 205-07° (decomp.) 206-07° (decomp.) 145°° 
Hydrochloride 206° (decomp.) “3 205-07° (decomp.) 
Sulphate 274° (decomp.) 274°8 (decomp.) 274° (decomp.) 
Picrate 234° (decomp.) 232°8 (decomp.) 232° (decomp.) 
Acetone compound 168° 168°20 
Product from treatment 284° (decomp.) 265°8 284° (decomp.) 
with dimethyl sulphate 
Tetrahydrobase 171° 213-15 171°! 
Dihydro base 164° 8Charred without 164°2 
melting 
TABLE II 
The red base Neprotine Jatrorrhizine 
(present work) (Chatterjee’s work) m.p. 
m.p. m.p. 
Base .. 206° (decomp.) Charring' without melt- 
ing 
Hydrochloride 205-08° (decomp.) Charrs! above 200° 206° (decomp.) 
without melting 

Nitrate 225° (decomp.) Be 225°4 (decomp.) 
Picrate 210-211° (decomp.)  208-10°* (decomp.) 
Tettahydrobase ~ 206°° 
0-Ethyltetrahydro base 124° 120°25 
O-Methyltetrahydro base 146° 146716 
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EXPERIMENTAL 
Plant material was supplied by Messrs. Ghosh and Co., Darjeeling. 


Isolation of the alkaloids of Mahonia nepalensis D.C.—Air-dried powdered 
roots and stem bark of the plants (2-5 kg.) were treated with alcohol (5-1). 
After two days, the yellow solution was drawn off and distilled on a water- 
bath. The brown residue was treated with water (500 c.c.), filtered, and the 
filtrate was treated with hydrochloric acid (congo red) and allowed to stand 
overnight. Crystals of the yellow base hydrochloride separated and were 
filtered and well washed, until the washings were yellow in colour. The 
brown filtrate was concentrated on a water-bath to 400 c.c. and filtered twice. 
On leaving this solution overnight, a second crop of yellow crystals appeared. 
The solution was decanted carefully into a shallow dish, evaporated partially 
and was left to crystallise. After two days red granules appeared and the 
dish was left undisturbed for a week till the quantity of granules stopped 
increasing. The red granules were filtered, washed and dried. Further 


spontaneous evaporation of the filtrate gave only a sticky brown mass which 
could not be crystallised. 


Yellow base hydrochloride—The foregcing yellow hydrochloride was 
crystallised first from hot water containing a few drops of hydrochloric acid. 
Recrystallisation from alcohol gave the hydrochloride (4-8 g./kg. of plant 


material) as yellow silky needles, m.p. and mixed m.p. with berberine hydro- 
chloride, 206° (decomp.). 


Yellow base-—The purified hydrochloride (4 2.) was dissolved in hot 
water and a dilute solution of sodium carbonate was added until alkaline, 
whereupon the yellow colour of the solution changed to dark brown, and 
on cooling, yellow silky needles of the free alkaloid (3 g.) were obtained. 
Repeated crystallisation from alcohol gave the base, m.p. 205-07° (decomp.). 
Found: C, 59-0; H, 5-9; N, 3°4. C.9H,;s0;N.3 H,O requires C, 59-0; 
H, 6:4; N, 3-4%. The melting point was undepressed by admixture with 
the higher melting form of berberine prepared as below :— 


Berberine sulphate (B.D.H.) was suspended in water and basified by 
adding saturated sodium carbonate solution. The solid was filtered and 
crystallised once from alcohol when yellow crystals of berberine, m.p. 145° 


were obtained. A second crystallisation from alcohol changed the m.p. to 
205-07°. 


The sulphate was prepared as usual and crystallised from alcohol had 
m,p, and mixed m.p. with terterine sulphate (B.D.H.) 274° (decomp.), 
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Picrate-—The picrate was prepared as usual and was crystallised from 
alcohol. It had m.p. and mixed m.p. with berberine picrate, 234° (decomp.). 


The acetone compound.—The yellow base hydrochloride (0-2 g.) dissolved 
in water (15 c.c.) was treated with sodium hydroxide solution (10%; 4 c.c.) 
and heated to 50°. Acetone (10c.c.) was added and the solution set aside 
when lemon-yellow needles of the acetone compound separated, which 
melted at 168°. Mixed with berberine-acetone compound, there was no 
depression in the melting point. 


Treatment of the yellow base with dimethyl sulphate-——An excess of di- 
methyl sulphate was added drop by drop to a solution of the yellow base 
(0-1 g.) in 10% potassium hydroxide solution with constant shaking. After 
half hour, yellow, silky fine needles separated. The shaking and the addi- 
tion of dimethyl sulphate were continued and the mixture was kept alkaline. 
After one hour, the crystals were filtered, well washed with water, dried and 
then recrystallised from alcohol as silky yellow needles, m.p. 284°, un- 
depressed by admixture, with the product, m.p. 284°, obtained from berberine 
by a similar procedure (Found: C, 55-0; H, 4-2%). 


Dihydrobase.—To a suspension of lithium aluminium hydride (0-8 g.) 
in ether (100 c.c.), the yellow base hydrochloride (0-5 g.) was added during 
the course of 15 minutes, and the stirring was continued for 12 hours. The 
complex was decomposed as usual and repeatedly extracted with ether. After 
the removal of ether, the dihydrobase was obtained as yellow crystals 
(0:25 g.). After recrystallisation from ether, the dihydro compound had m.p. 
and mixed m.p. with dihydroberberine, 164° (Found: C, 70°8; H, 5+3. 
CosH;g9O,N requires C, 71:2; H, 5-6%). 


Tetrahydrobase.—A solution of the yellow base hydrochloride (0-5 g.) 
in absolute alcohol (40 c.c.) containing platinum oxide (0-2 g.) was hydro- 
genated at 50 lb./sq. in. After the removal of solvent, the hydrochloride was 
obtained as a white powder, which crystallised from alcohol, had m.p. 
215-17°. The free base was obtained after basification with liquor ammonia. 
After recrystallisation from alcohol, the tetrahydro compound was obtained 
as white needles, m.p. and mixed m.p. with tetrahydroberberine, 172° 
(Found: C, 70-6; H, 6-2. CzoH,,0,N requires C, 70-8; H, 6-2%). 


Oxidation of the yellow base with potassium permanganate-—The base 
hydrochloride (0-5 g.) was dissolved in warm water (20 c.c.) and treated with 
potassium carbonate (0-15 g.) and then oxidised by slowly running in a solu- 
tion of potassium permanganate (2-5 g. in 50c.c. of water). The decolour- 
isation of the permanganate was very rapid at first and then slowed down 
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appreciably after 30 c.c. of the solution had been added. The contents were 
warmed on a water-bath for 3 hours, and left aside overnight. The mixture 
was filtered and the residual manganese dioxide was’ extracted twice with hot 
water and the combined filtrate concentrated to 30c.c. The solution was 
acidified with concentrated hydrochloric acid (congo red) and extracted con- 
tinuously with ether for several hours. After the removal of solvent, the 
gummy residue left behind was sublimed in vacuo. The sublimate was dis- 
solved in ethanol (5c.c.) and anhydrous ethyl amine (3 c.c.) was added and 
the solution was evaporated on a water-bath. The residue was heated at 
180° for 5 minutes and then sublimed in vacuo. The sublimate was crystal- 
lised from dilute alcohol to yield N-ethylhemipinimide as white, slender 
needles, m.p. and mixed m.p. with an authentic specimen, 90°. 


The red base hydrochloride—The red hydrochloride from the extraction 
of Mahonia nepalensis D.C. was crystallised from hot water containing a 
few drops of hydrochloric acid; the crystals were filtered and recrystallised 
from alcohol. Orange red prisms, m.p. 206-08° (decomp.), were obtained on 
slow crystallisation (0-2 g./kg. of plant material) (Found: C, 58-8; H, 
5:5. CopH220;NCI-H,O requires C, 58:6; H, 5-9%). 


The red base.—The purified hydrochloride (0-8 g.) was dissolved in hot 
water and a dilute solution of sodium carbonate was added until alkaline, 
when the brown colour of the solution changed to blood red and on coc ling, 
garnet red needles of the free alkaloid were obtained. It was recrystallised 
from alcohol to give the red base (0-6 g.), m.p. 206° (decomp.) (Found: C, 
59-1; H, 6-1. C.spH.,0;N-3H,O requires C, 58-7; H, 6:6%). The nitrate 
prepared by the usual method and crystallised from alcohol had m.p. 225° 


(decomp.). The picrate prepared by the usual method was crystallised from 
alcohol andhad m.p. 210-11°. 


Oxidation of the red base with potassium permanganate.—The red base 
(0-2 g.) on oxidation and working up in the usual way yielded hemipinic 
acid characterised as N-ethylhemipinimide, m.p. 90°. 


Tetrahydro base—A solution of the red base hydrochloride (0-5 g.) in abso- 
lute alcohol (30c.c.) containing platinum oxide (0-2 g.) was hydrogenated 
at 50 1lb./sq. in. After the removal of solvent, tetrahydrobase hydrochloride 
was obtained as a yellow powder which crystallised from alcohol as yellow 
needles, m.p. 206° (Found: C, 61:0; H, 6-7. C.ssH,,O0,NCI.H,O requires 
C, 60-7; H, 6°5%). 


The free tetrahydrobase was obtained by basification of the above hydro- 
chloride with liquor ammonia and was crystallised from alcohol, as white 
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needles, m.p. 206-07° (Found: C, 70-7; H, 6:7; N, 4-1. CooH,,0,.N 
requires C, 70:5: H, 6:7; N, 4:1%). 

O-Methyltetrahydrobase.—To a solution of the tetrahydrobase above 
(50 mg.) in absolute methanol (15 c.c.), a solution of diazomethane in ether 
(50c.c.) (prepared from 5g. of nitrosomethylurea) was added. Next. day, 
the solvent was evaporated and the residue was crystallised from absolute 
ether to give tetrahydropalmatine, m.p. and mixed m.p. with an authentic 
specimen, 146° (Found: C, 70:8; H, 5-7. C,,H,,0,N requires C, 70-9; 
H, 6:0%). 


O-Ethyl tetrahydrobase-—The tetrahydrobase (50 mg.) was ethylated 
in the same way as above using diazoethane. The product was crystallised 
from absolute ether, to give the ethyl ether as white needles, m.p. and mixed 
m.p. with d/-O-ethyltetrahydrojatrorrhizine (see below); 124°. 


dl-O-Ethyltetrahydrojatrorrhizine.—Il-Corypalmine (50 mg.) was ethy- 
lated in the same way as above using diazoethane. The product was dissolved 
in chloroform-methanol, treated with iodine in excess, and heated on a steam- 
bath with a little sodium acetate. The solvent and excess iodine were then 
removed and the residue treated on a steam-bath with zinc and dilute hydro- 
chloric acid until all of it had dissolved to a colourless solution. An excess 
of ammonia was added and the liberated base extracted with ether and the 
ether extract evaporated. The residue was crystallised repeatedly from 
ether to give colourless needles melting at 124°. 


SUMMARY 


Chemical investigation of the alkaloids of Mahonia nepalensis D.C. 
has established the identity of ‘ umbellatine ’ with berberine and ‘ neprotine’ 
with jatrorrhizine. 
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